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ABSTRACT 

Giant low surface brightness (GLSB) galaxies are commonly thought to be massive, dark matter dominated systems. However, this 
conclusion is based on highly uncertain rotation curves. We present here a new study of two prototypical GLSB galaxies: Malin 1 and 
NGC 7589. We re-analysed existing H I observations and derived new rotation curves, which were used to investigate the distributions 
of luminous and dark matter in these galaxies. 

In contrast to previous findings, the rotation curves of both galaxies show a steep rise in the central parts, typical of high surface 
brightness (HSB) systems. Mass decompositions with a dark matter halo show that baryons may dominate the dynamics of the inner 
regions. Indeed, a "maximum disk" fit gives stellar mass-to-light ratios in the range of values typically found for HSB galaxies. 
These results, together with other recent studies, suggest that GLSB galaxies are systems with a double structure: an inner HSB 
early-type spiral galaxy and an outer extended LSB disk. 

We also tested the predictions of MOND: the rotation curve of NGC 7589 is reproduced well, whereas Malin 1 represents a challenging 
test for the theory. 

Key words, dark matter - Galaxies: structure - Galaxies: kinematics and dynamics - Galaxies: individual: Malin 1 - Galaxies: 
individual: NGC 7589 - Methods: data analysis 



1. Introduction 

Low surface brightness (LSB) galaxies span a wide range of 
galaxy sizes, masses and morphologies, from dwarf spheroidals 
and dwarf irregulars to medium-size late-type disks (Sc-Sd) 
and to bulge-dominated early-type spirals (Sa-Sb). According to 
Beijersbergen et al. (1999), LSB spirals form a LSB Hubble se- 
quence, parallel to the classical HSB one. 

Giant low surface brightness (GLSB) galaxies are usu- 
ally considered as extreme cases of early-type LSB spirals 
(Beijersbergen et al. 1999). They have exceedingly extended 
LSB disks, with scale lengths ranging from ~10 to ~50 kpc, and 
central luminous components resembling the bulges of HSB spi- 
rals (e.g. Bothun et al. 1987; Sprayberry et al. 1995). They are 
massive systems, with L ~ L*. Like the ordinary LSB spirals, 
GLSB galaxies have low H I surface densities (Pickering et al. 
1997), but they are among the most gas rich galaxies known, 
with M m * 1O 1O M (Matthews et al. 2001). 

Few dynamical studies exist on GLSB galaxies. Pickering 
et al. (1997) studied four GLSB galaxies and found slowly ri- 
sing rotation curves, similar to those of late-type LSB disks (e.g. 
de Blok & McGaugh 1997; Kuzio de Naray et al. 2006). They 
concluded that GLSB galaxies are "the first examples of galaxies 
that are both massive and dark matter dominated". However, 
they warned that their rotation curves are highly uncertain, due 
to the low signal-to-noise ratio and the low spatial resolution of 
the observations. 

A slowly rising rotation curve in the presence of an inner 
concentration of light, as seen in those galaxies (see top panels of 
figure 7), is in marked contrast with the rule that there is a close 
correlation between the concentration of light and the shape of 
the rotation curve (Sancisi 2004 and references therein). For this 



reason, Sancisi & Fraternali (2007) started a re-analysis of the 
21-cm data from Pickering et al. (1997). This paper concludes 
that preliminary work with the study of the structure and dyna- 
mics of two GLSB galaxies: Malin 1 and NGC 7589. 

2. Hi Data Analysis 

We re-analyzed the H I datacubes of Malin 1 and NGC 7589, 
obtained with the VLA by Pickering et al. (1997). The angu- 
lar resolutions are low: 20". 9 x 20". 8 for Malin 1 (correspond- 
ing to ~ 32 kpc for an angular-diameter distance Da = 322 
Mpc) and 21".4 x 17".3 for NGC 7589 (~ 12 x 10 kpc for 
Da = 123 Mpc). We used the Groningen imaging processing 
system (GIPSY) (van der Hulst et al. 1992). The main physical 
properties of both galaxies are listed in table 1 . 

Total H I maps were obtained by adding the channel maps 
with line emission (from 24581 km s _1 to 24955 km s _1 for 
Malin 1 and from 8566 km s" 1 to 9157 km s" 1 for NGC 7589). 
In order to define the areas of emission in each channel map, 
we used masks, which were obtained by smoothing the original 
datacubes to lower spatial resolutions. We inspected by eye each 
channel map overlaying different masks. We obtained the best 
results by smoothing to 40" and 50" and clipping at 2cr s and 
3<x s for Malin 1 and NGC 7589 respectively (<x s is the rms noise 
in the smoothed cube). The resulting total HI maps are shown 
in figure 1 . 

From the masked data, we calculated the H I fluxes (correct- 
ing for primary beam attenuation). For Malin 1, we obtained 
S to t = 2.0 + 0.2 Jy km s , consistent within the errors with 
the value derived by Pickering et al. (1997) (S tot = 2.5 + 0.2 
Jy km s ) and the single-dish measurement by Matthews et al. 
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Table 1 : Physical Properties of Malin 1 and NGC 7589. Table 2: Kinematical Parameters of Malin 1 and NGC 7589 



Property 



Malin 1 



NGC 7589 



Ref. 



Parameter 



Malin 1 



NGC 7589 



Morphological Type 

AGN activity 

M B (mag) 

Redshift 

D L (Mpc) 

D A (Mpc) 

HI fluxdykms- 1 ) 

M H , (10 10 M o ) 

SFR FUV (M yr- 1 ) 

SFR NUV (M G yi-') 



SBO/a 
LINER 
-22.1 ±0.4 
0.0826 ± 0.0017 
377 ± 8 
322 ±6 
2.0 + 0.2 
6.7 + 1.0 
1.2 
2.5 



SAB(rs)a 1, 2 

BLAGN 1, 2 
-20.4 ± 0.2 3 

0.0298 ± 0.0017 4 

130 ±8 4 

123 ±7 4 

3.8 ±0.3 4 

1.5 ±0.3 4 

0.7 5 

1.1 5 



References. (1) Barth (2007); (2) Nasa Extra-galactic Database (NED); 
(3) Pickering et al. ( 1 997) give M v = -22.9 ± 0.4 for Malin 1 and Mr = 
-21.9 ± 0.2 for NGC 7589; Impey & Bothun (1989) give B - V = 0.8 
for Malin 1 ; Galaz et al. (2006) give B - R = 1 . 5 for NGC 75 89 ; (4) This 
work; (5) Boissier et al. (2008). 



(2001) (S, ot = 1.80 + 0.50 Jy km s" 1 ). For NGC 7589, we ob- 
tained Stat = 3.8+0.3 Jy km s _1 , whereas Pickering et al. (1997) 
derived 5 t ot = 2.7 + 0.3 Jy km s -1 . The differences are proba- 
bly due to the different masks used; the values are consistent 
within 2cr. Single-dish observations gave intermediate values, 
consistent with our value within about lcr: S tot = 2.96 + 0.51 
Jykms- 1 (Matthews et al. 2001)andS tot = 3.07+0.36 Jy km s _1 
(Springob et al. 2005) (this latter value is without the correction 
for HI self absorption). HI masses were estimated using the 
luminosity distance Z\, given by: 



D l = 7T b + o.5x(i- go )r]. 

"0 



(1) 



We assumed Hq = 70 km sec -1 Mpc -1 , £2 m = 0.27 and Qa = 
0.73 (qo = -0.59). Since we are dealing with low redshift ga- 
laxies, the exact value of go affects our results only slightly. The 
redshift was calculated as z — V sys /c, where V sys is the galaxy 
systemic velocity derived from the tilted ring analysis (see sec- 
tion 2.2). Our results are listed in table 1. 

2.1. Hi velocity fields 

The derivation of the velocity fields is not straightforward be- 
cause of beam-smearing effects, which are particularly severe 
here because of the low spatial resolution. Indeed, especially in 
the inner regions, the 21 -cm line profiles are unusually broad 
and asymmetric. A Gaussian fit or an intensity-weighted mean 
(IWM), as used by Pickering et al. (1997), give velocities far 
from the profile peaks and biased towards the systemic velo- 
city. This leads to velocity fields with systematically underes- 
timated inner velocity gradients and consequently to slowly ri- 
sing rotation curves. This is illustrated by figure 7 of Pickering 
et al. (1997), which shows the rotation curves of Malin 1 and 
of NGC 7589 overlaid on slices of the datacubes along the ma- 
jor axes of the galaxies. Clearly, the inner points of the rota- 
tion curves are at lower velocities than the emission peaks, sug- 
gesting that they are underestimated. The importance of beam- 
smearing effects is demonstrated by the galaxy models presented 
in section 2.4. 

We have derived new velocity fields with a different tech- 
nique to minimize beam-smearing effects. We selected veloci- 
ties near the peaks of the emission, estimating an IWM velocity 
from the upper part of the line profile (above ~ 75% of the peak 



a op ,(J2000.0) 
<5 opt (J2000.0) 
a km (J2000.0) 
4i„(J2000.0) 
V sys (kms" 1 ) 
Position Angle 
Inclination Angle 



12"36'"59 ! .4: 
14°19'49".3 : 
12 A 36 m 59\2 : 
14°19'51".3 : 



24766.7 ± 4.0 
see figure 3 

38° ± 3° 



23"18'"15 1 .7 ± P. 8 
00°15'40".2± 1".8 
23''18'"15 I .5 ±0 J .2 
00°15'37".3±2".5 
8926.4 ± 3.8 
307°. 7 ± l°.l 
58° ±3° 



Notes. The optical centres (a ovt , d opt ) are taken from NED. The inclina- 
tion angles are taken from Moore & Parker (2006) for Malin 1 and from 
Pickering et al. (1997) for NGC 7589. 



intensity) and neglecting the broad wings. This is like taking the 
profile peak velocities, but with the advantage of less noise in 
the resulting velocity field. In the outer regions of the galaxies, 
where the line profiles are more regular and symmetric, we used 
a Gaussian fit. The resulting velocity fields are shown in figure 1 . 

The left panels of figure 1 show the R-band image of Malin 1 
from Moore & Parker (2006) (top), the total H I map (middle) 
and the HI velocity field (bottom). The HI distribution is ap- 
proximately as extended as the optical LSB disk. The outer H I 
radius measures ~ 64", which corresponds to ~ 100 kpc. On the 
southern side, a spiral arm is visible. It extends for more than 80 
kpc and apparently it is the gaseous counterpart of the spiral arm 
detected by Moore & Parker (2006). The velocity field shows 
a strong variation of the position angle of the major axis with 
radius, both in the approaching and in the receding sides of the 
galaxy. Such behaviour is usually attributed to the presence of a 
warp. Overall, the velocity field is regular and symmetric. 

The right panels of figure 1 show the R-band image of 
NGC 7589 from Galaz et al. (2006) (top), the total H I map (mid- 
dle) and the H I velocity field (bottom). The H I distribution is as 
extended as the optical LSB disk. The outer H I radius measures 
~ 80", which corresponds to ~ 48 kpc. The HI distribution fol- 
lows that of the light, with two main concentrations (North-East 
and South- West) and two faint outer spiral arms. The velocity 
field is regular and symmetric. 

2.2. Rotation curves 

Kinematical parameters and rotation curves were derived by fit- 
ting a tilted-ring model to the observed velocity fields (Begeman 
1987). We assumed purely circular motions and described the 
galaxy by a set of concentric rings. We used all the points of the 
velocity fields to maximize the statistics of the least-square fit. 
The points were weighted by cos 2 (ff), where 9 is the azimuthal 
angle in the plane of the galaxy. 

As a first step, we determined the kinematical centre 
(akin, #kin) and the mean systemic velocity V sys . Then, we kept 
(akin,#kin) and V sys fixed and determined the position angle 
P.A. as a function of radius. Subsequently, we tried to deter- 
mine the inclination angle i. Unfortunately, because of the low 
spatial resolutions, inclination angles and rotational velocities 
are strongly correlated and no independent determination of ; 
was possible. Therefore, we kept the inclination angles con- 
stant, adopting values derived from optical observations. We 
used i = 38° for Malin 1 (Moore & Parker 2006) and i = 58° 
for NGC 7589 (Pickering et al. 1997). Inclinations derived from 
the total HI maps are consistent with those values. Finally, we 



2 



Lelli et al,; Structure and Dynamics of Giant Low Surface Brightness Galaxies 



1ALIN 1 




12"37 m 4 



2" CT 36"'5E 

RA (J2000.0) 



56 b 



NGC 7589 



23 h 1 8 m 1 5 s 



R band 



- U 15 



h 



Hl map 



- 17 




HI velocity field 




1 7 



0°15 



17 



- 0°15 



23 n 18 m 15 s 

RA (J2000.0) 



11* 



Fig. 1: Left Panels: Malin 1. Top: R-band image (Moore & Parker 2006). Middle: total HI map. Contours range from 0.9 x 10 20 to 
5.4 x 10 20 atoms cirT 2 with steps of 0.45 x 10 20 atoms ctrT 2 . Bottom: HI velocity field. The receding side is in dark-grey. Contours 
range from 24623 km s _1 to 24887 km s _1 with steps of 24 km s _1 (2 times the channel spacing). Right Panels: NGC 7589. Top: 
R-band image (Galaz et al. 2006). Middle: total HI map. Contours range from 1.1 x 10 20 to 9.9 x 10 20 atoms cm with steps of 
1.1 x 10 20 atoms crrT 2 . Bottom: HI velocity field. The receding side is in dark-grey. Contours range from 8772 km s _I to 9080 
km s with steps of 22 km s (the channel spacing). 3 
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Fig. 3: Left Panel: Position- Velocity diagram for Malin 1 along the Line Of Nodes (LON). The white squares show the projected 
rotation curve. Full lines show iso-emission contours from l<x to 8<x with steps of lcr, where lcr = 0.24 mJy/beam. Dashed lines 
show contours at -lcr and -2cr. Top-Right Panel: Position angle of the major axis for Malin 1 as a function of radius. The dotted 
and dash-dotted lines show the fits to the tilted-ring results. Bottom-Right Panel: Malin 1 Velocity Field. Iso-velocity contours range 
from 24647 km s _1 to 24887 km s _1 with steps of 24 km s _1 (2 times the channel spacing). The black line shows the LON. 
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Fig. 2: H I rotation curves for Malin 1 (top) and NGC 7589 (bot- 
tom), derived separately for the receding (up-triangles) and ap- 
proaching (down-triangles) halves. 



determined the circular velocity V c in each ring, keeping all the 
other parameters fixed. Our results are listed in table 2 and 3. 

In the tilted-ring analysis of NGC 7589, we used a ring width 
of 20" ~ 1 beam. For the kinematical centre, we obtained values 
which are consistent within the errors with the optical ones. The 
values obtained for V sys and P.A. approximately agree with those 
of Pickering et al. (1997). 

For Malin 1, we used a ring width of 10" ~ 1/2 beam. 
To determine the kinematical centre, we used the points inside 
40" ~ 62 kpc. The derived values agree within the errors with 
the optical ones. The value of V sys agrees within the errors with 
that derived by Pickering et al. (1997). The P.A. for Malin 1 is 
not constant. For the inner 15" ~ 24 kpc we found -1°.0 + 1°.6 
and for the outer radii a linear increase up to 65°. 1 + 4°. 9 at 
R = 65" (see the right panels of figure 3). For comparison, 
Pickering et al. (1997) let the P.A. vary from 0° + 10° to 34° + 3°. 
To derive the rotation curve of Malin 1 , we used a ring width of 
16" ~ 2/3 beam. In this way, the points of the rotation curve are 
almost independent. 

The errors a>. on the circular velocities have been estimated 



as o% = <j* 



+ CT 



asym 



, where is the formal error given 
by the fit and cr asym is an additional uncertainty due to asym- 
metries between approaching and receding halves. Following 
Swaters (1999), we define <x asym = (V Cjapp - V c ,rec)/4, where 
Vc,app an d Vc irec are respectively the approaching and receding 
rotation curves. Thus, we are assuming that the differences bet- 
ween the rotation curve from the entire velocity field and those 
for the approaching or receding side correspond to a deviation of 
2cr asym . Other uncertainties come from possible variations of the 
inclination angle. For example, for a change of i equal to 3°, the 
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Fig. 4: H I rotation curves for Malin 1 (left) and NGC 7589 (right). Diamonds show the rotation curves derived by Pickering et al. 
(1997), filled circles show our new rotation curves. The rotation curve of Malin 1 from Pickering et al. (1997) has been rescaled to 
; = 38°, adopted in our analysis. The open circle was derived only from the approaching half of Malin 1. 



rotation curve varies by about 14 km s 1 for Malin 1 and 5 km 
s- 1 for NGC 7589. 

Figure 2 compares the rotation curves derived separately for 
the approaching and receding sides. For NGC 7589, there is a 
high degree of symmetry between the two halves. For Malin 1, 
the symmetry is good out to R < 48" ~ 75 kpc. At R - 56" ~ 87 
kpc there is a difference of ~ 45 km s _1 between the two halves, 
whereas the last point at R = 72" ~ 112 kpc has been derived 
only for the approaching half. It has a large uncertainty and will 
be neglected. 

Figure 3 shows the line of nodes of Malin 1 overlaid onto 
the velocity field (bottom right panel) and the correspond- 
ing position-velocity diagram extracted from the datacube (left 
panel). The position- velocity diagram is not completely symme- 
tric: the approaching side is more extended than the receding one 
(see also figure 2). The projected rotation curve (white squares) 
closely follows the gas distribution. In particular, in contrast with 
figure 7 of Pickering et al. (1997), our first point of the rotation 
curve corresponds to the peak of the emission. 



2.3. Comparison between the new and the old rotation 
curves 

Figure 4 compares our new rotation curves with those derived by 
Pickering et al. (1997). The main difference is that our rotation 
curves rise more steeply and remain approximately flat out to the 
last measured points. As pointed out in section 2.1, the slow rises 
of the old rotation curves are due to the strong beam smearing 
effects in the inner regions. In the next section, we construct 3- 
dimensional galaxy models based on the rotation curves from 
Pickering et al. (1997) and those derived here and demonstrate 
that the latter are the correct ones. 

In the rotation curve of NGC 7589, our last two points lie 
below those of Pickering et al. (1997). The kinematical models 
described in section 2.4 show that Pickering et al. (1997) over- 
estimated the rotation curve at large radii, but the reason is not 
clear. 



Table 3: Rotation Curves for Malin 1 and NGC 7589. 



MALIN 1 



NGC 7589 



Radius 


V c 


0> c 


Radius 


V c 




(") (kpc) 


(km/s) 


(km/s) 


(") 


(kpc) 


(km/s) 


(km/s) 


8.0 12.5 


236.9 


9.4 


10.0 


5.9 


200.9 


9.3 


24.0 37.5 


216.2 


2.6 


30.0 


17.7 


194.6 


3.7 


40.0 62.5 


197.4 


4.7 


50.0 


29.5 


204.8 


3.2 


56.0 87.5 


214.1 


11.8 


70.0 


41.3 


200.1 


3.3 



Notes. Columns (1) and (2): ring radius in arcsec and kpc. Column (3): 
circular velocity. Column (4): error on the circular velocity as defined 
in the text. 



2.4. Kinematical models 

In order to verify the correctness of the new rotation curves, 3- 
dimensional galaxy models were built. We used the GIPSY task 
GALMOD, that creates model datacubes by describing a galaxy 
disk with a set of rings. Each ring is characterized by: the kine- 
matical parameters (xo,yo), V sys , PA. and i, the rotation velocity 
V c , the gas velocity dispersion cr H i , the face-on gas column den- 
sity Zhi and the scale-height h z . 

We created two different sets of models. The first set was 
built out of our new kinematical parameters and rotation curves. 
These are the "steeply rising curve" models. The second set was 
built out of the kinematical parameters and rotation curves de- 
rived by Pickering et al. (1997). These are the "slowly rising 
curve" models. The HI surface density profiles were derived 
from the total H I maps, azimuthally averaging over a set of el- 
lipses, determined by the assumed kinematical parameters. For 
the H I vertical distribution, we assumed an exponential profile 
exp(-z/h z ), with h z = 200 pc. For the gas velocity dispersion, 
we used <x HI = 8 km s over the entire disk. The models are 
almost indistinguishable for differences in o"hi of ±4 km s . 
All the models were smoothed to the spatial resolution of the 
observations, thus they also reproduce beam-smearing effects. 
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Fig. 5: Position- Velocity diagrams extracted from the observations and the model datacubes. From left to right: observations; model 
using the steeply rising rotation curve and the other parameters derived in this work; model using the slowly rising rotation curve 
and the parameters derived by Pickering et al. (1997). For Malin 1 {top) we took a slice along P.A. = 35°, for NGC 7589 {bottom) 
along P.A. = 305°. Full lines show iso-emission contours from 2cr to 9<x with steps of lcr. Dashed lines show contours at -2cr. For 
Malin 1 {top) <x = 0.24 mjy/beam, for NGC 7589 {bottom) cr = 0.35 mJy/beam. 



In figure 5, we compare position-velocity diagrams from the 
models and the data. For NGC 7589, we took a slice through the 
datacubes along a position angle of 305°. This is a mean value 
between our result of 307°. 7 and that of 302° by Pickering et al. 
(1997). For Malin 1, we took a slice along a position angle of 
35°. This is a mean value for the major axis. 

For both galaxies, the observed datacubes (left panels of fi- 
gure 5) show the presence of H I emission at high rotational ve- 
locities near the galaxy centre {R < 20"). The "slowly rising 
curve" models (right panels) do not reproduce such emission, as 
the H I is spread from low rotational velocities near the centre 
to high rotational velocities at large radii. On the contrary, in the 
"steeply rising curve" models (middle panels) the H I emission 
is concentrated at high rotational velocities, as seen in the data. 
Both models show tails of emission toward the systemic velocity 
due to beam-smearing effects. Overall, the observed data are re- 



produced better by the "steeply rising curve" models based on 
our new results than by the "slowly rising curve" models based 
on those by Pickering et al. (1997). 

We also note that the "slowly rising curve" model of 
NGC 7589 at large radii (R > 50") exhibits HI emission at 
rotational velocities higher than those observed. This demon- 
strates that the rotation curve by Pickering et al. (1997) is over- 
estimated in the outer regions. 



3. Mass Models 

The steeply rising rotation curves found for Malin 1 and 
NGC 7589 suggest that GLSB galaxies have a dynamical beha- 
viour more similar to a HSB than to a LSB galaxy. To determine 
the relative contributions of luminous (gas and stars) and dark 



6 



Lelli et al,; Structure and Dynamics of Giant Low Surface Brightness Galaxies 



MALIN 1 

Radius (arcsec) 
20 40 




40 60 
Radius (kpc) 



3 « 



NGC 7589 

Radius (arcsec) 
40 60 




20 30 
Radius (kpc) 



3 « 



Fig. 6: HI surface density profile for Malin 1 (left) and NGC 7589 (right), corrected for inclination and cosmological dimming 



matter to the gravitational potential, we built mass models fol- 
lowing Begeman (1987). 

3.1. Gas and stars 

The contribution of the gaseous disk was computed using the 
surface density profiles derived from the total HI maps (fi- 
gure 6). These were multiplied by a factor (1 + z) 4 to correct for 
cosmological dimming and by a factor 1 .4 to take into account 
the presence of Helium. Molecular and ionized gas were not ex- 
plicitly considered in the mass model. However, since they are 
usually distributed as the stellar component, their contribution 
is reflected in a small increase of M*/L. Consistently with the 
models built in chapter 2.4, we assumed an exponential vertical 
distribution with a scale height of 200 pc. 

To compute the contribution of the stellar component, we 
used surface brightness profiles from the literature. For the 
vertical distribution of the stellar disk, we assumed Z(z) = 
sech 2 (z/zo)/zQ (van der Kruit & Searle 1981a,b), with zo = 
300 pc. 

For Malin 1, we used the I-band profile from Barth (2007) 
for R < 10 kpc and the R-band profile from Moore & Parker 
(2006) for R > 10 kpc. We estimated the R - 1 colour of the in- 
ner galaxy regions using the SDSS data and the relations from 
Fukugita et al. (1996). We found r' - i' = 0.46, which corre- 
sponds to R - I = 0.70. Using this colour, the match between the 
two profiles is very good (see the upper-left panel of figure 7). As 
noted by Barth (2007), Malin 1 is a "normal" SBO/a galaxy, em- 
bedded in an extended LSB disk. This double structure is clearly 
shown by our combined luminosity profile. The HSB disk ex- 
tends out to about 20 kpc, where a photometrically distinct LSB 
disk appears. 

For NGC 7589, we used the R-band profile by Galaz et al. 
(2006) for R < 25 kpc and the R-band profile by Pickering et al. 
(1997) for R > 25 kpc. In the inner regions, where both profiles 
are available, there is a systematic difference, maybe due to the 
different calibration of the observations. Therefore, we applied 
a correction of 0.2 mag to the points given by Pickering et al. 
(1997), by requiring that the two profiles match in the inner re- 
gions. This correction is within the Icr uncertainties in the pho- 
tometric measurements. The result is shown in the upper-right 
panel of figure 7. Inside ~ 20 kpc, the luminosity profile is typi- 



cal of a HSB galaxy with a bulge (R < 5 kpc), a lens component 
(R ~ 5 - 13 kpc) and a HSB exponential disk (R ~ 13 - 20 kpc). 
In the outer regions, a photometrically distinct LSB disk extends 
out to ~ 50 kpc. 

3.2. Dark matter 

For the dark matter distribution, we assumed pseudo-isothermal 
haloes described by equation: 



where po is the central density and rc is the core radius, po and 
rc are free parameters of the mass models. 

We also tried rotation curve decompositions with NFW dark 
matter haloes (Navarro et al. 1995) but found that the results are 
not significantly different. 

3.3. Rotation curve decompositions 

The H I rotation curves have much lower spatial resolution than 
the surface brightness profiles (and their resulting dynamical 
contributions). Thus, a least-square fit done only on the observed 
points of the rotation curves would neglect some features present 
in the luminosity profiles. To avoid this problem, we have li- 
nearly interpolated the rotation curves and performed a fit with 
steps of 1". Errors on the rotation curves were used as weights 
in the fit. 

Figure 7 shows the "maximum disk" decompositions of the 
rotation curves. The resulting mass-to-light ratios are M»/Lr = 
3.4 forMalin 1 and MJL R = 2.5 for NGC 7589 (Table 4). These 
values are somewhat uncertain due to the lack of resolution in 
the central regions. However, we estimate that the mass-to-light 
ratio cannot be higher than 5 for Malin 1 and 3 for NGC 7589. 
Thus, they surely are in the range of values typically found for 
HSB early-type spirals (Noordermeer 2005). Under the maxi- 
mum disk hypothesis, the baryons dominate the galaxy dyna- 
mics out to a radius of R ~ 20 kpc for Malin 1 and R ~ 30 kpc 
for NGC 7589. In contrast with the results of Pickering et al. 
(1997), dark matter does not necessarily dominate everywhere. 

Both Malin 1 and NGC 7589 show colour gradients, which 
may indicate different stellar populations. Bothun et al. (1987) 



7 



Lelli et al,; Structure and Dynamics of Giant Low Surface Brightness Galaxies 



MALIN 1 

Radius (arcsec 
20 40 



NGC 7589 

Radius (arcsec) 
40 60 




40 60 
Radius (kpc) 



20 30 
Radius (kpc) 



Fig. 7: Upper Panels: R-band surface brightness profiles for Malin 1 (left) and NGC 7589 (right). See text for details. Lower Panels: 
rotation curve decompositions for Malin 1 (left) and NGC 7589 (right). Dots show the observed rotation curves. The curves show 
the contributions due to gas, stars and dark matter and the resulting total rotation curve. 



Table 4: Mass Models for Malin 1 and NGC 7589. 



Parameter 


MALIN 1 


NGC 7589 


M,/L R (M Q /L G ) 
Po (10 -3 M Q pc- 3 ) 
rc (kpc) 

XI- 


3.4 

121.4 ±93.8 
2.2 ±0.9 
4.6 


2.5 
2.3 ±0.3 
19.9 ± 2.0 

4.0 



Mhsb/^r (M Q /L e ) 3.3 2.6 

M L sb/Lr (M q /L q ) 0.9 0.9 

p (10 -3 M G pc -3 ) 73.9 ±36.5 2.4 ± 0.2 

r c (kpc) 2.9 ±0.8 20.7 ± 1.6 



Notes. Top: mass models using a single stellar component and maximi- 
zing the stellar mass-to-light ratio. Bottom: mass models using a HSB- 
LSB decomposition. The stellar mass-to-light ratio of the LSB disk was 
fixed at 0.9, whereas that of the HSB disk was maximized. See text for 
details. 



studied the luminous "blobs" in the outer disk of Malin 1 (see 
figure 1) and found a mean B - V colour of 0.37, whereas in the 
inner parts (R < 20") the mean B - V colour is 0.90 (Impey & 
Bothun 1989). Galaz et al. (2006) traced the B - R colour profile 
of NGC 7589 out to R ~ 20 kpc and found complex trends, but in 
general the galaxy tends to become bluer at large radii. It makes 



sense, therefore, to decompose the luminosity profiles and use 
two different M s /Lr for the inner and the outer regions. 

We separated the LSB disk by fitting an exponential profile 
to the points beyond R > 20 kpc. The resulting disk parame- 
ters are: pu(Q) = 24.7 mag arcsec -2 , h = 51.7 kpc for Malin 1 
and p R (0) = 23.3 mag arcsec -2 , h = 14 kpc for NGC 7589. 
The HSB component was isolated by subtracting the LSB disk 
from the total photometric profile (see top panels of Fig. 8). We 
maximized the contribution of the HSB component and found 
Mh S b/L r = 3.3 for Malin 1 and M hsb /£r = 2.6 for NGC 7589, 
very similar to the previous values. Maximizing also the contri- 
bution of the LSB disk, the outer parts of the rotation curve can 
be reproduced without the need for a dark matter halo. However, 
the resulting Mlsb/^r ratios would be unrealistically high (~25 
for Mahn 1 and -15 for NGC 7589), as it is usually found by 
maximizing the disks of "ordinary" LSB galaxies (e.g. de Blok 
et al. 2001). Therefore, we estimated Mlsb/^r = 0.9 using the 
relations of Bell et al. (2003) and assuming B - V = 0.44. The 
latter value is the mean colour found by McGaugh & Bothun 
(1994) for a sample of bulgeless LSB disks. The results of the 
fits are listed in table 4 and shown in the bottom panels of fi- 
gure 8. 

We also performed a standard bulge-disk decomposition for 
the inner HSB part and found a degeneracy of the mass-to-light 
ratios of bulge and disk. Different combinations of Mb u i ge /iR 
and Mdisk/^R provide comparable fits, as the rotation curves are 
not resolved in the inner regions. By maximizing both bulge and 
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MALIN 1 NGC 7589 
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Fig. 8: Upper Panels: R-band surface brightness profiles for Malin 1 {left) and NGC 7589 (right). The curves show the HSB-LSB 
decompositions. See text for details. Lower Panels: rotation curve decompositions for Malin 1 (left) and NGC 7589 (right). Dots 
show the observed rotation curves. The curves show the contributions from the HSB and the LSB stellar disks, the gaseous disk, the 
dark matter halo and the resulting total rotation curve. 



disk we found that their contributions to the first point of the ro- 
tation curve are approximately equal. Therefore, in this "maxi- 
mum light" hypothesis, the inner steep rise is not explained only 
with the bulge, but also the disk is dynamically important. This 
further suggests that Malin 1 and NGC 7589 have an inner "nor- 
mal" HSB disk (see also Barth 2007). 

Following McGaugh (2005), we used the derived values of 
M t /L R and M m to check the position of Malin 1 and NGC 7589 
on the baryonic Tully-Fisher relation. We found that both ga- 
laxies follow the relation within the observed scatter. 

3.4. MOND 

The modified newtonian dynamics (MOND) was proposed by 
Milgrom (1983) as an alternative to dark matter. This theory 
modifies the force law at accelerations a lower than a critical 
value flo and naturally explains the flatness of rotation curves 
(see the review by Sanders & McGaugh 2002). 

GLSB galaxies are ideal systems to test MOND. Firstly, they 
are very diffuse and extended, thus the test is possible down 
to very low values of a/ao (deep MONDian regime). Secondly, 
they have massive gaseous disks. In the MONDian framework, 
this gives an important contribution to the galaxy dynamics, 
which is not subject to uncertainties as that of stars (i.e. the value 
of M*/L). For the MOND fit to the rotation curves, we followed 
Begeman et al. (1991) and assumed Oq = 3000 km 2 s~ 2 kpc 1 



(Bottema et al. 2002). Two different M/L ratios for the inner and 
outer regions of the galaxies have been used as free parameters 
in the fit. 

For NGC 7589, we obtained a good fit by using a bulge- 
disk decomposition. The stellar mass-to-light ratios predicted 
by MOND are acceptable: M buige /L = 4.7 and M disk /L = 1.3. 
For Malin 1, it is preferable to use the HSB-LSB decomposition 
shown in figure 8, but it is not possible to leave both Mhsb/L and 
Mlsb/L free. Indeed, the fit would require a negative value for 
Mlsb/Lr. This suggests that the presence of the extended LSB 
disk, combined with the large amount of gas, creates problems 
for MOND. We fixed M LSB /L R and used only M H sb /Lr as a 
free parameter. Figure 9 shows the result for Mlsb/^r = 0.5. 
Using the relations by Bell et al. (2003), this corresponds to 
an extremely blue colour for the outer disk, i.e. B-V = 0.2. 
Higher values of Mlsb/Lr worsen the fit, whereas lower va- 
lues only slightly improve it, as the gas becomes dominating in 
the outer regions. A bulge-disk decomposition, as that used for 
NGC 7589, does not improve the fit. Higher values of the critical 
acceleration, as ao = 3700 km 2 s~ 2 kpc -1 (Begeman et al. 1991) 
or a = 4000 km 2 s~ 2 kpc -1 (McGaugh 2004), worsen the fit. 
Interpolation functions different from the standard one, as those 
proposed by Famaey & Binney (2005) or Milgrom & Sanders 
(2008), give slightly worse fits. 

Although the rotation curve of Malin 1 is not of high qua- 
lity, the discrepancies between the observed velocities and the 
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MALIN 1 - MOND FIT 




Radius (kpc) 

Fig. 9: MOND applied to Malin 1 using two different M t /L R for 
the inner HSB region (Mhsb/^r = 3.7) and the outer LSB disk 
(Mls,b/Lr = 0.5). Dots show the observed rotation curve. Lines 
show the Newtonian contributions due to the gaseous, the HSB 
and the LSB disk and the result of the MOND fit. See text for 
details. 



MOND prediction are significant as they are equal to ~ 25 km 
s _1 , i.e. ~ 12 — 13%. However, the warp in the outer disk of 
Malin 1 may greatly affect these results. The rotation curve used 
here was derived by varying the P.A. with radius (see figure 3) 
but assuming a constant inclination angle of 38°. It was not pos- 
sible to obtain a better determination of i from the velocity field. 
A change of the inclination angle in the outer parts of only 6° 
from i = 38° to i = 32°, which cannot be ruled out, would result 
in a rotation curve in total agreement with the MOND predic- 
tion. Note that here we have used the simplest formulation of 
MOND (Milgrom 1983), that may not strictly apply to warped 
non-axisymmetric disks. 

4. Discussion 

4.1. GLSB galaxies and their double HSB-LSB nature 

It is known that there is a dichotomy between HSB and LSB 
systems (e.g. Verheijen & Tully 1999). HSB galaxies usually 
have steeply rising rotation curves and can be described by a 
maximum disk, whereas LSB galaxies have slowly rising ro- 
tation curves and are thought to be dominated by dark-matter 
everywhere. The dynamics of Malin 1 and NGC 7589 is typi- 
cal of HSB systems. The concentration of mass, indicated by 
the steeply rising rotation curve, follows the inner concentration 
of light. For both galaxies, maximum-disk solutions have been 
found with reasonable values of M*/X. Thus, we conclude that in 
the inner regions either baryons are dominating or, in case dark 
matter dominates, it must follow the distribution of light. This 
may be in line with some recent numerical simulations (Tissera 
et al. 2009), which suggest a gravitational effect of baryons on 
shaping the distribution of dark matter. 

Both surface photometry and gas dynamics indicate that 
Malin 1 and NGC 7589 have a double structure: an inner "nor- 
mal" spiral galaxy (with a bulge and an HSB disk) and an outer 
extended LSB disk. Is such a double structure present in the 
other GLSB galaxies as well? The other two GLSB galaxies 
studied by Pickering et al. (1997) are F568-6 (Malin 2) and 
UGC 6614. F568-6 is kinematically lopsided and difficult to 
interpret. UGC 6614 was observed in Ha by McGaugh et al. 



(2001) who found a steeply rising rotation curve (see their fi- 
gure 5) and thus a dynamical behaviour similar to that of a HSB 
system (see also de Blok et al. 2001). Walsh et al. (1997) stu- 
died the GLSB galaxy NGC 289 using optical and radio observa- 
tions. They found an inner HSB disk and a steeply rising rotation 
curve that can be fitted by a maximum disk with M*/L\ = 2.3. 
They drew conclusions very similar to ours. Another example is 
NGC 5383. This barred spiral galaxy fits in the GLSB class. It 
has a LSB disk with h ~ 9.7 kpc and a G-band absolute magni- 
tude of M G = 20.6 (van der Kruit & Bosma 1978) 1 . Its surface 
brightness profile is very similar to that of NGC 7589 and ex- 
tends out to ~ 54 kpc (see Barton & Thompson 1997). The HI 
mass is Mui = 5.9 x 10 9 M and the gas dynamics is typical of 
a HSB system (Sancisi et al. 1979). 

Sprayberry et al. (1995) catalogued another 13 objects as 
GLSB galaxies. As argued by Barth (2007), their inner HSB 
disk could have been missed due to the low spatial resolution 
of the optical observations. However, for some of them, HI 
rotation curves were derived. These galaxies are NGC 5533 
(Noordermeer et al. 2007), NGC 5905 (van Moorsel 1982), 
NGC 4017 (van Moorsel 1983), NGC 2770 (Garcfa-Ruiz et al. 
2002), UGC 2936 (Pickering et al. 1999) and PGC 45080 (Das 
et al. 2007). All of them show a dynamical behaviour typical 
of HSB systems, with the possible exception of PGC 45080 
that is kinematically lopsided (see figure 1 1 of Das et al. 2007). 
Therefore, it seems unlikely that GLSB galaxies constitute a sub- 
class of the bulge-dominated LSB galaxies (e.g. Beijersbergen 
et al. 1999, Galaz et al. 2006). It is likely, instead, that they are 
"normal" HSB galaxies which possess an outer extended LSB 
disk. 

GLSB galaxies do show some physical properties closer to 
HSB than to LSB systems. Usually, LSB galaxies are quite 
amorphous and rarely host an AGN (Impey & Bothun 1997; 
Bothun et al. 1997). In contrast, GLSB galaxies can have 
a well-defined spiral pattern (Beijersbergen et al. 1999) and 
show nuclear activity with the same probability as HSB spi- 
rals (Schombert 1998). Moreover, LSB galaxies are deficient in 
molecular gas (e.g. O'Neil et al. 2003; Matthews et al. 2005), 
whereas CO emission has been detected in some GLSB galaxies 
(e.g. O'Neil & Schinnerer 2004; Das et al. 2006). 

At this point it is natural to ask how common it is for HSB ga- 
laxies to have diffuse outer LSB stellar disks. This question has 
been addressed in the past (see e.g. Bosma & Freeman 1993 and 
references therein) and examples were found such as NGC 3642 
(Verdes-Montenegro et al. 2002). More recently, Erwin et al. 
(2008) studied a sample of 66 SO-Sb galaxies and found that 
~ 24% of them show "anti-truncated" disks, i.e. the surface 
brightness profile of the disk can be described by two exponen- 
tials, the outer one with larger scalelength than the inner one. 
Such "up-bending" luminosity profiles have also been observed 
in late-type spiral galaxies (Hunter & Elmegreen 2006; Pohlen 
& Trujillo 2006). GLSB galaxies may be extreme examples of 
anti-truncated early-type HSB spirals. Moreover, GALEX obser- 
vations have recently revealed extended XUV-disks around HSB 
spirals (Thilker et al. 2007). Also some GLSB galaxies have 
been observed with GALEX (Boissier et al. 2008) and, at such 
wavelengths, they resemble the XUV-disk galaxies of Thilker 
et al. (2007) (see in particular Malin 1, Malin 2 and NGC 7589). 

Regarding gaseous disks, it is known that a large number of 
spiral galaxies, perhaps the majority, have HI disks extending 
far beyond the bright optical ones (e.g. Sancisi et al. 2008). The 



1 The original values have been rescaled to a distance of 33.7 Mpc, 
as we assume Hq = 70 km s~' Mpc -1 . 
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outlying HI has typical surface densities of 1-2 M e pc~ 2 , close 
to the values measured in GLSB galaxies. Similar results have 
been found also for early-type galaxies in the field. Oosterloo 
et al. (2007) have observed at 21 -cm a sample of 30 gas rich 
E/SO galaxies and discovered that about 2/3 are embedded in ro- 
tating H I disks with masses, densities and sizes comparable to 
the gaseous disk of Malin 1 . Deep optical observations would 
be useful to verify whether such spirals and ellipticals have ex- 
tended LSB stellar disks associated with the gaseous ones. If this 
is the case, these objects could be considered to be in the same 
class of GLSB galaxies. 



4.2. The formation of GLSB galaxies 

GLSB galaxies are a challenge for theories of galaxy formation. 
The main properties to be explained are: i) the double HSB-LSB 
nature, ii) the large extent of the outer LSB disks and iii) the 
regular and symmetric large-scale kinematics of the latter, as ob- 
served in both Malin 1 and NGC 7589 (see figure 2). Given the 
very large orbital periods at the outermost radii (e.g. fd yn ~ 2.8 
Gyr for Malin 1 and t dyn ~ 1.4 Gyr for NGC 7589) it is clear 
that these outer disks must have been in place and undisturbed 
for several Gyrs. 

CDM cosmological simulations of galaxy formation tend to 
produce dense and compact disks, as much angular momentum 
is lost during the hierarchical assemblage (see Kaufmann et al. 
2007 and references therein). Hoffman et al. (1992) showed that 
rare density peaks (~ 3cr) in underdense environments (voids) 
may lead to the formation of a GLSB galaxy. However, GLSB 
galaxies show the same clustering properties of HSB spirals 
and are not necessarily associated with voids (Sprayberry et al. 
1995). 

Mapelli et al. (2008) simulated the evolution of a collisional 
ring galaxy. The propagation of the ring leads to the formation 
of a large and diffuse system, with structural properties similar 
to those of GLSB galaxies. However, these simulations predict 
slowly rising rotation curves and therefore a different dynamics 
from that of GLSB galaxies. 

Noguchi (2001) proposed a secular-evolution scenario, start- 
ing from a massive barred HSB galaxy. The bar-instability can 
redistribute the matter and the angular momentum in the disk, 
but the scale length h cannot increase by more than a factor of 2- 
2.5 (see also Debattista et al. 2006). Therefore, this mechanism 
cannot explain the formation of an extended LSB disk as that of 
Malin 1 (with /; ~ 50 kpc). Moreover, it should explain why the 
majority of barred early-type spirals do not have an outer LSB 
disk (Erwin et al. 2008). 

Simulations by Penarrubia et al. (2006) investigated the for- 
mation of LSB disks around HSB galaxies by interactions with 
dwarf companions. They found that the dwarf galaxies are dis- 
rupted by tidal forces and the debris settle on an extended disk 
with a nearly exponential profile. Interestingly, the scale length 
of the resulting disk depends on the compactness of the original 
stellar distribution in the dwarfs and, in extreme cases, it can be 
as large as ~ 50 kpc. Thus, it would be possible to explain a wide 
range of disk structural parameters as those observed in GLSB 
galaxies. However, such simulations predict a decrease in the 
rotation velocity of the outer stellar disk of about 30-50 km s 
with respect to the inner one. Such a decrease is not observed in 
the HI rotation curves of GLSB galaxies. 



5. Conclusions 

We studied the structural and dynamical properties of two gi- 
ant low surface brightness galaxies: Malin 1 and NGC 7589. For 
both galaxies, we re-analyzed existing 21 -cm line observations, 
obtained with the VLA by Pickering et al. (1997), and derived 
new rotation curves. The quality of these curves was tested by 
building model datacubes. Finally, we decomposed the rotation 
curves using available optical surface brightness profiles and as- 
suming pseudo-isothermal dark matter haloes. We also studied 
the predictions of MOND. Our main results can be summarized 
as follows: 

1. in contrast with the previous results by Pickering et al. 
(1997), both Malin 1 and NGC 7589 have steeply rising ro- 
tation curves. From a dynamical point of view, these GLSB 
galaxies are more similar to early-type HSB than to late-type 
LSB disks. This agrees with the result of the optical study 
by Barth (2007) that described Malin 1 as a normal SBO/a 
galaxy surrounded by a huge LSB disk. Both surface pho- 
tometry and gas dynamics indicate that also NGC 7589 has 
such a double HSB-LSB structure. 

2. In these two GLSB galaxies, the concentration of mass in- 
dicated by the steeply rising rotation curve corresponds to 
the inner concentration of light. The rotation curves of both 
galaxies can be fitted by a maximum disk with M*/Lr ~ 3. 
Contrary to the case of the less massive LSB disks where 
dark matter is believed to dominate everywhere, in the inner 
regions of GLSB galaxies either the baryons dominate or the 
dark matter follows the distribution of light. 

3. Malin 1 provides a severe test for MOND because of its ex- 
traordinary size and high H I mass. The rotation curve de- 
rived here differs significantly from the one predicted by 
MOND. However, the discrepancies would disappear if the 
outer disk of Malin 1 were warped and turned to more face- 
on. Unfortunately, the present observations are not adequate 
to determine the trend of the inclination angle. 
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